Antiferromagnetism, Stripes, and Superconductivity in the 
t-J Model with Coulomb Interaction 
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We study mean-field phases of the t-J model with long-range Coulomb interaction. In the order 
of increasing doping density we find a classical antiferromagnet, charge and spin stripes, and a 
uniform d-wave superconductor, at the realistic doping parameters. Both in-phase and anti-phase 
stripes exist as metastable configurations, but the in-phase stripes have a slightly lower energy. The 
dependence of the stripe width and the inter-stripe spacing on the doping is examined. Effects of 
fiuctuations around the mean-field states are discussed. 



The cuprate materials which exhibit high-Tc supercon- 
ductivity show many different ordering tendencies as the 
hole doping concentration {x) in the system is varied. At 
zero doping the cuprates are antiferromagnetic insulators 
below the Neel temperature (Fig. 1(a)) The ordered 
antiferromagnetic phase ceases to exist with more than 
about 2% of holes. Around 5% doping the system be- 
gins to show superconductivity at low temperatures. In 
the intervening 2 — 5% the system is lacking either anti- 
ferromagnetic order or superconductivity. In the doping 
range between 5% and 15%, superconductivity co-exists 
with the (dynamic) stripe order [H-|4|. The stripes are 
most readily seen in neutron scattering experiments at 
non-zero energy transfer, hence "dynamic" , where one 
finds evidence of one-dimensional periodic modulation of 
the antiferromagnetic order as well as the charge density 

i- 

As noted from the early days of high-Tc, the cuprates 
are characterized by strong electron-electron repulsion 
1^, which gives rise to a Mott-insulating state at half- 
filling. The t-J model, which incorporates such electron 
repulson, has been extensively studied in connection with 
high-Tc superconductivity. Although it is quite likely 
that the t-J model correctly captures the short-distance 
correlation of electrons in the cuprates, one must not a 
priori overlook the long-range part of Coulomb interac- 
tion. In particular for the t-J model without Coulomb 
interaction, phase separation occurs for a wide range of 
doping concentration , which pre-empts the possibility 
of a high-pairing-scale superconducting state. 

Various mean-field theories which assumes a uniform 
ground state have been in existence These theo- 

ries have varying degrees of success in understanding 
the phases of the cuprates. Such mean-field approaches 
have however been subject to the skepticism that the 
no-double-occupancy constraint is treated only approxi- 
mately. Attempts to improve the treatment of the con- 
straint result in a strongly fluctuating gauge fleld. Re- 
cently one of us (D.-H.L.) was able to integrate out the 
gauge field exactly at low energies in the uniform su- 
perconducting phase It was shown that despite the 
drastic modification of the excitations, mean-field vacua 



serves as a good starting point in characterizing the zero- 
temperature state. 

This work is performed under the assumption that 
mean- field theory will capture the short-distance/high- 
energy ordering tendency of the t-J -I- Coulomb model, 
and that the long-distance/low-energy properties can be 
understood by studying the soft fluctuations of the mean- 
field order parameters. 

The Hamiltonian we consider is the following: 



H = -tY^{b,b\flho.+ h.c.) + JY^{S. ■ S,- -n,n,) 

{ij) {ij) 



2 ^ r, 

i^j 



The no-double-occupancy constraint is expressed as 
fiafia+blbi — 1 = in terms of the spinon (fia) and holon 
(bi) operators. Other notations include = 1 — blbi (site 
electron density), and fi = 1 — x (average electron den- 
sity). Our first goal is to understand the mean-field 
phases sustained by this model. Unlike other mean- 
field theories in the past we include the magnetic 
order parameter on an equal footing with all other or- 
der parameters. Among other things this gives us the 
possibility of obtaining the long-range ordered antiferro- 
magnetic state observed in experiments. 

Mean-field theory 

Our mean-field theory, in essence, is a variational ap- 
proach: a trial wavefunction is constructed and parame- 
ters are varied to obtain the minimum energy. We con- 
struct a wavefunction which allows local magnetic mo- 
ments (but not limited to antiferromagnetism), super- 
conducting pairing (but not limited to d-wave symme- 
try), and modulations in the charge density. The trial 
wavefunction j^E') is given by j^*) — \'^b) ® 1^*/) where 
the bosonic (l^'h)) and fermionic (l^*/)) states are in- 
dependently constructed from their respective vacua as 
follows: 
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I*/) = n ("."/It + ^i/.i) (^?/u + ^?/'t) |0/). (2) 

a,b 

Repeated indices j and I implies summation over lat- 
tice sites. The bosons are assumed to be condensed, 
and the fermion ground state is constructed by occu- 
pying the (yet undetermined) quasiparticle orbitals la- 
beled by a,b. The mean- field single-particle wavefunc- 
tions Uj,Vj, w^, z\ and Xj ^-re varied to minimize 

(vj/iifivf) - + ^lJ^c. - m 

i 

-^Y,l^^\n,-n\^). (3) 

i 

Lagrange multipliers Ai, and /i are introduced to guar- 



antee that the average occupation obeys the constraints 
locally as well as globally. 

The calculation is carried out numerically on a x Ny 
lattice with a periodic boundary condition. Not assum- 
ing any translational invariance, we first perform totally 
unrestricted minimization for N^^Ny < 16. After the 
nature of the solution is established, we perform a more 
restricted search for Nx,Ny up to 120. The results re- 
ported below are for t/J — 3, Vc/J — 5. Other choices of 
t/J and Vc/J are also studied, with results that are not 
qualitatively different from those presented below. 

We find three prominent types of order. In the order of 
increasing doping, they are antiferromagnetic insulator, 
charge/spin stripes, and uniform d-wave superconductor. 
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FIG. 1. Mean-field phases of the t-J-|-Coulomb model: (a) antiferromagnet (b) in-phase stripe (c) anti-phase stripe (d) 
uniform d-wave superconductor with residual antiferromagnetic order (e) non-magnetic uniform d-wave superconductor. Yel- 
low/red circle represents sites fully /partially occupied by the electrons. Magnetic moments for partially occupied sites (white 
arrow) are smaller than those of fully occupied sites (black arrow). The bonds shown in black and green are the pairing 
amplitudes, which differ in sign for x and y directions. The dashed vertical lines in Fig. 1(c) indicate the boundary where the 
shift of the antiferromagnetic order parameter occurs. 



Antiferromagnet (Fig. 1(a))- At zero doping the 
mean-field ground state shows antiferromagnetic long- 
range order. Each electron is surrounded by four neigh- 
bors with opposite spins. It is an insulator because the 
occupation constraint forbids the electrons to hop. This 
mean-field state is a caricature of the insulating anti- 
ferromagnet observed in the undoped cuprates |^]. For 
x<Xci w0.02, the doped holes are localized, often in the 
form of elongated puddles (or finite-length stripes). How- 
ever, these localized puddles do not disrupt the overall 
antiferromagnetic order. 
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FIG. 2. The stripe spacing Ac vs. doping for the in-phase 
stripes. The numbers above the arrow indicate the approxi- 
mate doping fraction where the periodicity changes. 



Stripes (Fig. 1(b)- (c)): For Xd < x < Xc ^ 0.14 the 
mean-field ground state shows charge corrugation in the 
form of stripes. A stripe is a region which is extended 
in one direction (say y) and localized in the other, with 
partially occupied sites. There are two types of stripes, 
anti-phase and in-phase. The antiferromagnetic order 
parameter goes through a n phase shift across the anti- 
phase stripe, whereas it remains in-phase across the in- 
phase stripe. As a result the anti-phase stripe modulates 
the antiferromagnetic order with a period twice that for 
the charge density. 



We find that the in-phase stripes are the ground states 
for Xci < X < Xc- At X = Xc Su first order phase transi- 
tion occurs, after which the system becomes uniform. In 
all cases we have studied the in-phase stripes are bond- 
centered and have a width of two lattice constants. The 
stripe spacing, on the other hand, depends on doping 
and increases as the doping level decreases. The smallest 
distance between the in-phase stripes we observe is four 
lattice constants, and it occurs near Xc. This spacing is 
maintained in the doping range 0.095 < x < 0.14. For 
this range the site hole density in each stripe varies from 
0.19 per site (corresponding to a line density of 0.38) at 
X = 0.095 to 0.28 (line density 0.56) at a; = 0.14. As 
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the doping decreases below 0.095 a new stripe configu- 
ration with the stripe spacing equal to five lattice con- 
stants emerges as the ground state. The stripe width is 
still two. The evolution of the stripe spacing with the 
doping concentration is shown in Fig. 2. We infer from 
this the existence of a series of integer stripe spacings 
as doping decreases. Each spacing has a non-zero range 
of stability giving rise to plateaus in the modulation pe- 
riod of the charge density. The step-wise evolution of 
the stripe spacing is clearly a lattice commensuration ef- 
fect. In the presence of thermal fluctuation of the shape 
of stripes (quantum fluctuation is known not to roughen 
the stripe shape), we expect a smoother evolution of the 
modulation period. Inside each stripe there exists strong 
superconducting pairing as well as weak magnetism, as 
illustrated in Fig. 1(b). The pairing gap inside the in- 
phase stripe is comparable to the maximal pairing gap 
(« 0.05J) observed in the uniform c?-wave phase, while 
the magnetic moments are a fraction of the full moment, 
(Sz) = ±1/2, of the insulator. An in-phase stripe is in 
some sense an optimum pairing state which is confined 
in the one-dimensional geometry. As in-phase stripes 
get closely spaced, it is likely that transverse fluctua- 
tion smears out the charge corrugation and results in a 
high-pairing-scale superconductor. 

In the entire range of Xci < x < Xc the anti-phase 
stripes are metastable mean-field solutions. The energy 
difference between the most favorable anti-phase and in- 
phase stripes is shown as function of doping in Fig. 3. 
The largest difference (10%) occurs at x = 0.025 and the 
smallest (0.14%) at Xc- Note that anti-phase stripes come 
very close in energy to the in-phase ones near Xc, and 
therefore, fluctuations that are omitted by the mean-field 
theory may stabilize the anti-phase stripes. (One such 
candidate is the transverse fluctuation of the stripes.) 
When that happens, the progression of the ground states 
vs. doping will be according to Figs. 1(a) through 1(e). 
We discuss the properties of the anti-phase stripe in the 
following. 

The anti-phase stripes are also bond-centered and have 
a width equal to two lattice constants. The stripe spacing 
evolves in a step- wise fashion similar to Fig. 2. The range 
of hole density inside the anti-phase stripe is consequently 
similar to the in-phase stripe case. The anti-phase stripes 
also have non-zero pairing and magnetic moments (Fig. 
2(c)). However, the pairing scale is considerably smaller 
(by a factor of three) than that in the in-phase stripes. 
In this sense, anti-phase stripes are unfavorable as far as 
pairing is concerned. 

Uniform d-wave superconductor (Fig. l(d)-(e)): The 
homogeneous phase, x > Xc, is characterized by d-wave 
pairing and, for x near Xc, some residual antiferromag- 
netism JTo| . The pairing scale is maximum dX x — Xc 
where /S.max ~ 0.05J and decreases monotonically as 
X increases. The antiferromagnetic moments disappear 
completely for x > 0.2. 
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FIG. 3. Energy difference of the anti-phase and the 
in-phase stripes, 2{Ea„ti- Ein)/{Ea„ti + Ei„), as a function 
of X. In-phase stripe has an increasingly lower energy as the 
average hole density diminishes. 



We find it gratifying that our mean-field theory pro- 
duces states which extrapolate between extreme classical 
(antiferromagnet at a; = 0) and extreme quantum (uni- 
form superconducting) limits. 

Coulomb interaction and high pairing scale: In light 
of the above mean- field results, we argue that high-Tc 
superconductivity is a cooperative effect due to the short- 
range magnetic exchange and the long-range Coulomb 
interaction. [ pd| Larger antiferromagnetic exchange fa- 
vors higher pairing scale, however it also causes phase 
separation to set in at a lower doping level and the 
uniform high-pairing state becomes inaccesssible. At 
t/ J = 3, our model shows a phase-separated ground 
state for x < 0.26 in the absence of Coulomb interaction. 
Coulomb interaction suppresses phase separation and 
produces two compromises - stripes and a high-pairing- 
scale d^-2_yi superconductor. 

Fluctuations 

The low-energy excitations often appear in the form of 
the fluctuation of the order parameters. In the present 
context these include: phase fluctuation {9b) of the Bose 
condensate phase fluctuation {9p) in the pairing 

condensate of f), orientation fluctuation of the mag- 
netic moments (O), shape fluctuation and displacement 
of the stripes, gapless (neutral) fermion excitation in the 
case of d-wave pairing, and most importantly the "gauge 
fluctuation" inherent in the slave-boson theory ||l^ . 

Due to the occupancy constraint and the form of 
Eq. (0) there exists an internal gauge symmetry [p^ un- 



der a local phase change, hi 



piSiU. f. 



e'^'f^a- Such 



symmetry is broken by most of the mean-field vacua, giv- 
ing rise to a fluctuating gauge field as the soft mode. 
Recently one of us were able to integrate out the gauge 
fluctuations exactly in the non-magnetic, uniform d-wave 
superconducting phase corresponding to Fig. 1(e) 
The result is the confinement of two Goldstone modes, 
9(, and 9p, into one 4> = 29i, — 9p, which is the phase 
of the electron superconducting condensate. The fi- 
nal low-energy dynamics is that of a phase-fluctuating 
superconductor with gapless fermion excitations |13|,p[: 
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£ = £0 + -C^- + ^int, where 



2 

£^ = ^ -Ipnaidt - iVxnTxdx " iVynTzdy)lpn 
n=l 



(4) 



In the above p is the average Cooper pair density, 
is the fermion field associated with the nth gap 



node, Tr is the third PauU matrix, 



specify the 



hnear dispersion of the gapless fermions, and j^^ = 

\{Y.n^'n<yTzi^nc,,iVxli'la-^la,iVy2-lp2a'lp2a), IS the fcrmion 



3-current. Due to the gauge fluctuation the parameters 
u, are strongly renormalized. In particular K is 
proportional to x, which accounts for the low superfluid 
density in spite of the high pairing scale in the pseudogap 
regime. Similar treatment of gauge fluctuations has been 
done for each of the mean-field phases discussed above. 
The results are somewhat technical and will be reported 
elsewhere. 

For the antiferromagnet the mean-field state satisfies 
the occupation constraint exactly. In this phase the only 
low-energy degree of freedom is the fluctuation of the di- 
rection of the local spin. The spin degrees of freedom 
is gauge-neutral and hence unaffected by strong gauge 
fluctuations. The interaction between the spin waves is 
described by the familiar non-linear sigma model 
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(5) 



In two space dimensions, it is well known that the spin- 
wave fluctuation does not destroy the antiferromagnetic 
long-range order as long as the spin stiffness Ka- is suf- 
ficiently big. This certainly seems to be the case for the 
undoped cuprates 0]. 

In the uniform d-wave phase corresponding to Fig. 
1(d) there exists residual antiferromagnetic moments. 
The low-energy degrees of freedom are those of the non- 
magnetic superconductor plus the spin fluctuation. Since 
the wavevector associated with the commensurate anti- 
ferromagnetic ordering, k = (7r,7r), mismatches the mo- 
mentum connecting the gap nodes, the magnetic degrees 
of freedom decouple from the low-energy fermions. The 
effective action is then simply the sum of Eq. (||) and 
Eq. (^) . Due to the smallness of the magnetic moments in 
this phase, however, the quantum fluctuation is likely to 
wash out the long-range correlation of the residual mag- 
netism. In that case, the distinction between magnetic 
and non-magnetic superconductors becomes obscure. 

The fluctuations in the stripe phase is the richest. The 
low-energy degrees of freedom include the phase of the su- 
perconducting condensate and fermion quasiparticle exci- 
tations inside the stripes, the fluctuation of the magnetic 
moment, and the displacement and shape fluctuation of 
the stripes. Unlike the antiferromagnetic and the uniform 



superconducting phases, the fluctuation of stripes can not 
only modify the properties of a given stripe phase, but 
also change the energy ordering between the in-phase and 
anti-phase stripes. Despite some progress p^ , a com- 
plete picture which involves all of the above fluctuations 
is still lacking. The subject is currently under investiga- 
tion. 
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